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Abstract 

Motivated by a recent discovery of Supernova 2010gx and numerical results of Fryer et al. (2010), we 
simulate light curves for several type I supernova models, enshrouded by dense circumstellar shells, or "super- 
wind" , rich in carbon and oxygen and having no hydrogen. We demonstrate that the most luminous events 
like SN 2010gx can be explained by those models at moderate explosion energies ~ (2-^3) foe if the total 
mass of SN ejecta and a shell is ~ (3 -f- 5) Mq and the radius of the shell is ~ 10 16 cm. 

1 Introduction 

Recently, iFrver et all (l2010h have presented arguments in favour of type la supernova (SN) explosions taking 
place within rather dense extended C-0 envelopes formed during a merging process in Double Degenerate 
(DD) scenario. They obtained light curves which are very powerful in hard spectral range and last sometimes 
very long in visible lig ht. On the other hand, there are observations of a very luminous type Ic SN 2010gx 



(|Pastorello et a.l.ll2010f l . and it was suggested that they may shine due to a shock wave propagating in C-O- 



rich circumstellar material. 

Motivated by these results we construct light curves for several non-evolutionary type I supernova models 
which have dense C-0 shells or winds around them. 

2 Presupernova models 

We have not tried to compu te evolutionary or hydrodynamic models leading to formation of structures dis- 
cussed by Fryer et al. (|2010T). All our pres upernova models are constructed artificially as described elsewhere 



(jChugai et al.ll2004t iBaklanov et al.ll2005l) . Their main parameters are given in Table [TJ The interior part 



is a quasipolytropic model in mechanical equilibrium where temperature is related to density as T cx p 031 . 
This part has mass M e j and radius i? e j, the subscript "ej" denotes here that this part becomes a supernova 
ejecta after explosion. In case of type Ib/c SN simulations M e j can be much less than the total mass of the 
collapsing core and the condition of mechanical equilibrium is not necessary, it is just a convenient form of 
parameterization of models. The outer parts have power-law tails, p cx r _p , with p given in the Table [TJ Mm 
denotes the mass of radioactive 56 Ni in ejecta, and M w , i? w are mass and radius of the "wind". No attempt 
is done to keep equilibrium in the "wind" , but the dynamical times of those huge envelopes are so large that 
no appreciable motion has developed during the time of light curve simulation. 

The firs t set o f models tries to mimic what is described (though very tersely) in IFrver et~aL ( 2010l ). 



Fryer et all (|2010l ) claim that a double degenerate (DD) merger event leads to a density profile like p cx r 4 



which is shown in their Fig. 5, as well as a p cx r variant. However, this structure extends only up to 



3i?© there, while they see shock interaction up to r ~ 10 cm ~ 10 Rq in their Fig. 10 (jFrver et al 



2010). It is hard to imagine the formation of such an extended structure on a dynamical time-scale of DD 



event. Nevertheless, it is interesting to consider those extended structures ind ependently of DP sce nario 
because they can help to explain extremely powerful supernovae like SN 2010gx (IPastorello et alJl201Cft . We 



say a few words on other feasible ways to formation of these dense circumstellar shells in Sec. [5j 
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Fig. 1: Left: typical composition of "la" models used in simulations. "Fe" includes "Ni" (which is not only 56 Ni) and 
other iron peak elements. "C" and "O" lines almost coincide. Right: density profiles as a function of M r for steepla 
(p oc r~ 4 ; solid) and medium bigla (p oc r~ 3 ; dashed) models. 
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Fig. 2: Left: density profiles as a function of r for steepla (solid) and mediumla (dashed) models. (Model medium bigla 
has the same structure of the "wind" like the latter, p oc r~ 3 , but its outer radius is larger, R w — 3 ■ 10 6 i?©). Right: 
density profiles as a function of M r for the same models. 



We have two types of distribution of chemical composition in our simulations. A typical composition of 
the the first type is shown in Fig. [TJ It tries to mimic a distribution of elements which can be produced in a 
thermonuclear explosion like in SN la. The models with this composition have a suffix "la" in their names. 

The second type, "lb", has just a uniform distribution of elements in C-0 envelope like in outermost 
layers of "la" model (in left panel of Fig. [IJ. 

Density profiles for models with the "wind" p(r) oc r~ 4 and p(r) oc r~ 3 are shown in Fig. [T] All models, 
except for medium bigla, have T = 10 3 K in the "wind", and medium bigla has T = 2.5 • 10 3 K there. 
Already at T = 2.5 • 10 3 K we have got a spurious flash of light emitted by the huge envelope. 



3 Light curve simulations 



We perform calculations of the synthetic light curves u sing our multi-group radiation hydrodynamic code 



STELLA m its 



Blinnikov ct al 



stand ard 
20061) . 



setup (Blinnikov et all Il998l iBlinnikov fe Sorokinal 12004 iBaklanov et~"aT1 120051: 



The explosions have been simulated for all "la" variants by a "thermal bomb" with energy 1.6 TO ergs 
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Table 1: Models (all 


masses 


M and 


radii 


i? are 


in solar units) 


Model 


M cj 




M Ni 


P 


M w 


Rw 


E, foe 


steepla 


2 


1 


0.7 


4 


0.1 


3 • 10 5 


1.6 


medium bigla 


1.4 


1 


0.7 


3 


0.9 


3- 10 6 


1.6 


mediumla 


1.4 


1 


0.7 


3 


0.8 


3- 10 5 


1.6 


shallowlb 


1 


10 





2.5 


2.9 


10 5 


3 


standardlb 


0.2 


10 





2 


3.5 


8- 10 4 


3 


bright lb 


0.2 


10 





1.8 


4.8 


9- 10 4 


1, 2, 4 



1.6 foe like in iFrver et al.1 (|2010h . A similar "thermal bomb" was used for "lb" runs (the burst duration 
is 10 seconds in the innermost layers of ejecta, AM = 0.06Mq) with variable energy E, see Table [TJ In 
radiation hydrodynamics runs we used 200 (100 in the ejecta plus 100 in the "wind") radial mesh zones for 
"la" runs and 300 ones (150 plus 150) in "lb" runs. All runs employed 100 frequency group in transport 
solver and relatively short spectral line list (~ 1.5 • 10 5 ) in the opacity routine. 



4 Results 

We report the numerical results for a set of selected models below. 



4.1 Type la explosions within C— O shells 

The light curves for "la" runs are shown in Figs. [3H3 

The main effect which is clearly visible is the dependence of the flux and duration of the light curve on 
the outer radius in the mediumla (i.e. p{r) oc r~ 3 ) models. A physical explanation for this may be a longer 
diffusion time in the larger model. However, we should investigate also the dependence of the results on grid 
resolution. The presented "la" runs are done on grids with 200 radial by 100 frequency mesh points. We 
have checked them on cruder grids (90 radial by 100 frequency meshes). The behavior near maximum does 
not change much, while the tail part is more sensitiv e to the grid . 

A direct comparison of our results with those of Frver et all ( 2010h is not possible, because we do not 
know their exact initia l cond itions. Our mediumla models (p(r) oc r -3 envelopes) have about twice as much 



C-0 as in IFrver et al. I ll201Cft to exaggerate the effect of the cir cumstellar matter. Still for a smaller outer 



radius of the wind our light curves seem to evolve faster than in IFrver et al.1 (|2010f ). 

One should also note that in Fig. [4] we see initial "plateau" due to relatively high initial T — 2.5 • 10 3 K 
in the "wind" . The radiation energy stored there is simply emitted away. The process of formation of the 
extended structu re may be much lo nger than the time of its cooling. It seems that at least part of the early 
light in Fig. 9 of Frver et al. ( 2010l) may be attributed to this spurious emission. 



4.2 Type lb explosions within C— O shells 



Supernova SN 2010gx dPastorello et al. 2010) is extremely luminous and extremely interesting. For other 
very luminous events (jOfek et al.l 12007^ ISmith et al.l 120071 l20ld : iGal-Yam et~all 120091: lYoung et all l2010h 
some models involve explosions on a h ypernova scale, where a huge amount of 56 Ni is produced (see, e.g., 
Nomoto et aj]l2007t iMoriva et al.ll2010L for type Iln SN 2006gy and type Ic SN 2007bi). 

Our goal is to explain SN 2010g x on another way, with minimu m energy of explosion. This seems to be 
possible based on the old idea due to Grasberg fc Nadvozhin (1986): when we have two subsequent explosion 
events. The first explosion is weak and produces the dense circumstellar structure (which we call "wind" 
here, but which is not a steady wind, of course). The second, normal SN explosion, produces very bright 
light due to the shock embedded int o a dense medium. A physical mechanis m for those multiple explosions, 
proposed bv lHeger fc Woosievl ( 20021 ) (pulsation pair instability) was used bv lWooslev et al. ( 2007 ) to explain 
SN 2006gy with moderate energy of ~ 3 foe without any radioactive material. One should keep in mind that 
there may be other, still unexplored, routes to repeated explosions in stellar evolution, especially in binary 
systems. 
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Fig. 3: Light curves of the model steepla with p(r) oc r~ 4 in outer layers. 
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Fig. 4: Light curves of the model medium bigla with p(r) oc r~ 3 in outer layers with a larger outer cut of the wind 
(i? = 3-10 6 i? Q ). 
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It is interesting to look at the results of our "lb" runs with zero 56 Ni mass, where one can see a pure 
effect of radiative shock in C-0 envelope producing a very bright supernova. 

One clearly needs a very large radius to produce a bright event. The radius should not be too large, 
otherwise the diffusion of photons is too slow and the light curve becomes too wide. One needs high density 
for strong production of visible light by the shock, but not too high, otherwise the mass of the "wind" and 
the optical depth of the shell become too large. 

We have run tens of different models, varying all parameters. Here we show only a small part of our 
results, see Figs. EHZ1 Note, that the scale for the luminosity L has changed in comparison to Figs. [3] [5] The 
values of p — 4 and p — 3 in the p{r) oc r~ p distribution produce too steep density gradients for our goal 
and we try smaller values of p. The models denoted shallowlb have p — 2.5. We call models standardlb 
if they have p — 2 because this is a standard value for a steady stellar wind, not for our set of models. The 
best value we find is p = 1.8 in the brightlb set. 

There is an indirect evidence t hat the density pr ofiles around supernovae having strong circumstellar 



interaction may be s hallow indeed (iPrieto et al.l 120071) . or do have rather complicated structure due to pre- 
supernova evolution (jDwarkadas et al .I l2010h . 



4.3 Comparison with SN 2010gx. Role of opacity 

Our code Stella allows us to produce the light curve in different filter systems. He re we present them i n 
SDSS ug ri filters using the sam e standards and transmission functions as we did in I Phillips et al. I (120071) . 



Following iPastorello et al. (2010) we assume redshift z — 0.23 for SN 2010gx. However, we take the distance 



modulus=40. 28 which we find for this redshift in standard cosmology with Hq = 71 km/s/Mpc, Q m = 0.27, 
and J7a = 0.73. Pastorello et al. ( 201dl ) do not give the distance modulus explicitly. We can guess that they 



assume a larger distance, because they compare the values of observed and absolute g magnitudes for the 
host galaxy. The larger distance can be obtained in a flat universe only for very low Hq and 17 m , perhaps 
excluded by current observational data. 

Fig. [8] shows that our models brightlb easily reach observed maximum fluxes in g filters for various 
values of explosion energy. Energy E = 2 foe is perhaps too low, energy E = 4 foe is a bit too high, and the 
decay of light is too fast. The fluxes in uri filters are a bit lower than in observations, but before twiddling 
around those models we should note a very important factor, namely, line opacity. 

In our standard Stella setup we treat expansion opacity in lines as for type la supernovae, where we have 
homologous expansion and isotropic velocity gradient dv/dr — 1/i, with t time elapsed after the explosion. 
Now we have dv/dr 1/t, because the light is produced in the radiative shock attached to a very dense 
shell (see Sec. 14.41 below). To check the effect of line opacity we have run several tests when the expansion 
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Fig. 6: Light curves of the models shallowlb and standardlb. 
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Fig. 8: Observed and synthetic light curves in ugri filters of models brightlb with energies E = 2 foe (left) and E = 4 foe 
(right). 
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Fig. 9: Observed and synthetic light curves in ugri filters of models bright lb with energies E = 2 foe (left) and E = 4 foe 
(right) for higher expansion effect in line opacity. 



opacity is calculated with the fixed value dv/dr = 1 day . Fig. [9] shows that in this case the observed flux 
is higher in many bands. 



4.4 Hydrodynamic profiles of different models 

We present in Fig.[lU]a comparison of details in distributions of density, temperature, velocity, luminosity and 
Rosseland optical depth in models steepla and mediumla for the same date a few weeks after the explosion. 

We see that a shock wave in the steepla model is much stronger (due to a steeper density gradient) and 
leaves the grid quickly. The luminosity is due partly to heat release from the shocked matter and partly to 
56 Ni, but not to the shock itself: for the epoch shown in Fig. [10] the shock has already left the grid. Note 
that units of velocity on the plot are 10 9 cm/s, i.e. ten thousand km/s. Fo r a shallower profile p(r) oc r -3 the 



shock is "buried" deep and forms a dense shell like in our SN Iln models (jChugai et al.ll2004t IWooslev et al 



2003). One can see that the shock associated with this dense shell is contributing to lumin osity (although 



the 56 Ni role is not negligible). It is interesting to study models with zero mass of 56 Ni (as in lWoosley et al.l 
( 2007t n and we do this for SN 2010gx. 



We do not observe a box-like profile of temperature as in Fig.10 oflFW er et al. (2010) 



FigHU shows profiles for one of our models of SN 2010gx. We see a formation of very dens shell and 
luminosity production on the shock. 



5 Summary and discussions 



The aim of the current work was to check numerically the role of shock interaction in producin g light by shocks 
in ca rbon-oxygen "winds" around type I supernovae. T he problem arises from observations (jPastorello et al 



2010l and references therein), as well as from theory ( Frver et al. 201Clh . Our main conclusion is that the 
shock interaction is very important indeed and it cannot be ignored in many cases. 

We do not know exact initial conditions used bv lFrver et al. ( 2010h . and we do not confirm their conclusion 



on long-living shocks in steep (p(r) oc r ) density profiles. We suspect that the shock lives there due to 
another distribution of matter (there is perhaps a medium with p(r) = const on their grid). In part of our 
results we see initial flashes of light due to initial temperature in the "wind", even if it is as low as 2500 K. 
Nevertheless, we find long-living radiating shocks in all profiles with p(r) oc r~ 3 and in the less steep ones. 

A big question remains on the mechanism of the formation of those huge and dense envelopes. What is 
the time-scale of this process? How far can the envelope extend in radius? What is the density profile and 
the temperature of matter before the core explodes? 

If those structures form in reality then we can have extremely bright and long light curves. This question 
deserves further investigation. The effect of the extended envelopes is probably more important not for SNe la 
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Fig. 10: Profiles of density (black lines), velocity (in 10 9 cm/s, blue), luminosity (logarithm of absolute value in 10 40 ergs/s, 
violet), logarithm of matter temperature (green) and Rosseland optical depth (red). The scale for density is on the right 
Y axis, for all other quantities it is on the left Y axis. Left panel: model steepla with p(r) oc r~ 4 in outer layers. Right 
panel: model mediumla with p(r) cx r~ 3 in outer layers 
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but for very luminous SNe Ib/c discovered recently ( Pastorello et al.ll2010() . In that case one may think about 
pulsational pair instability as a means for formation of extended shell. One can also speculate about mergers 
of white dwarfs with CO-cores of WR stars. The merger event may lead to a moderate explosion with energy 
of a few percent of foe which forms a cloud of matter around the presupernova expanding with relatively low 
velocity. Then, if the core collapse follows within years after the merger, one can have all needed conditions 
for very bright supernovae. 

Our simulations show that the very luminous SN 2010gx may indeed be produced by a supernova explosion 
if there is enough surrounding material for a shock to transform the kinetic energy of ejecta into observed 
light. We find that SN 2010gx can be explained at explosion energies ~ (2 3) foe for a non-steep density 
profile, if the total mass of SN ejecta and a shell is ~ (3 -j- 5) M© and the radius of the shell is ~ 10 16 cm. 

The fits to fluxes in individual filters are not yet perfect in our simulations. This is natural: we have 
taken quite arbitrary and primitive chemical compositions, density distributions, etc. But in many cases the 
synthetic fluxes are higher than observations, and this looks encouraging. One can try building a better 
fit to observations by variations of initial conditions in the model. However, it seems that it is to early to 
optimize the models along these lines. This optimization will probably not give us a true insight and a better 
understanding of the problem. 

E.g., many technical subtleties still remain in the treatment of line opacity. First, we should correct for 
the expansion effect not only the flux equation, bu t also the energy equ ation. Microscopically, for exact 



monochromatic opacity, there is no expansion effect (IBlinnikov 19961 199711. For aver age absorption opacity 



there is some effect, but it can not be treated so simply as by Hoeflich et al. ( 19931 ). fsee lSorokina fc Blinnikov 
2002). Moreover, there is another complication with anisotropic velocity gradient. Thus, before optimizing 



the fits, one has to build new techniques for radiation transfer in these conditions. 

We have not discussed observed line spectra: for type Iln, by definition, one clearly sees narrow emission 



lines produced in the shells. Not so in SN 2010gx: narrow circumstcllar lines are not seen ([Pastorello et al 



2010h . There is no hydrogen, which is easily excited, and the most abundant elements, carbon and oxygen, 



should be present perhaps as C II and O II ions in the envelope (see the temperature in Figs. I11IT2"]) . These 
ions do no t have many strong lines in visible light. It is not easy to identify C and O lines on photospheric stage 
in SNe Ic ( Young et al.ll20l"of) . and now they should be excited even if there is no radioactive material. So one 
has to look for weak and narrow lines in noisy spectra. This problem certainly deserves further investigation 
with account of different conditions for ionization/excitation of shells under the shock radiation. 

The main complication to the whole picture is possible fragmentation of the dense shell. The attempts 
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on multi-D treatment of S N ejecta evolution are rather old (jTenorio-Tagle et alJll991 ; Chevalier fc Blondin 
199.4 iBlondin et aj]ll99fit). mo re recent results and references may be found in (|Dwarkadasl 120071 l2008h . 



See also Ivan Marie et ah ( 2010h for the case of SN 2006gy, but without real treatment of radiative transfer . 
There are several 3-D MC transport co des (|Hoeflichl 120021 iLucvl 120051: iKasen et aLll2006i 120071: ISiml 12007 : 
Tanaka et al. 20081 Kromer fc Siml 20091) but thev a re not actually coupled to hydrodynamics and there are 



many difficulties in doing this ( Almgren et al. 2010l ). 

Full NLTE treat ment is needed to predict spectra, but very little is done on this even for SN Iln. E.g., 
Dess art et al. (2009) are surprisingly successful in reproducing the spectra of SN 1994W in a set of atmo- 



spheric models, but their method is applicable only to monotonic velocity structures, not to shocked shells . 
Moreover, one should be cautioned about the relation of "photospheric" radius found bv lDessart et al. (2009) 
which shrinks, and the radius of the shocked shell in SNe Iln which grows. This is already explained by 
Smith et all (|2008L l2010h . 

Nevertheless, we conclude that, provided the formation of rather dense and extended circumstellar shells, 
the extremely powerful events of type Ib/c like SN 2010gx can be explained with moderate energy of explosions 
without invoking any radioactive material. 
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